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1. Introduction 

In high-energy physics, transverse-momentum-dependent parton distribution functions (TMD- 
PDFs), with or without spin dependence, have proved to be an essential quantities for unraveling 
the internal structure of protons |JIJ], as well as being an ingredients representing hadronic physics 
in a wide class of factorized physical observables. In this effort we will consider the unpolarized 
TMDPDF and how it should (or could) be defined in a way compatible with a generic factoriza- 
tion theorem of a transverse-momentum-dependent (TMD) physical observable. To properly define 
such quantities one needs to consider the following issues combined: 

1 . The role of the soft function, to be defined below, as a crucial part of the factorized hadronic 
tensor. 

2. The distinction between pure and naively calculated collinear matrix elements. 

3. How a definition of a TMDPDF should be independent of the use of any regulator of the 
non-ultra-violet divergences (nUV). By definition, "nUV divergences" include the physical 
infra-red (IR) ones of perturbative QCD (pQCD) and the un-physical ones like rapidity di- 
vergences. 

Consider the Drell-Yan heavy lepton pair production process where the lepton pair is produced 
with transverse momentum qj much larger than Aqcd- Assuming that the soft and collinear modes 
completely capture the IR of full QCD, which is a highly nontrivial statement, then, at leading 
twist, one can establish the following factorization theorem for g^-dependent cross section: 

o = H{Q 2 /y 2 )J n JnS , (1.1) 

where the g^-dependence is implicit in the pure collinear matrix elements J nJ! and the soft function 
S. By "pure" collinear contribution we mean that the soft (zero-bin) contamination has to be sub- 
tracted out [|h. In the remainder of this section we refer to all quantities in Eq. ( |L~I| ) as the partonic 
version of the physical matrix elements unless otherwise stated. To explain in simple terms how 
the TMDPDF should be defined we emphasize the following: 

1 . In pQCD the extraction of the hard part is performed by considering a generic Feynman 
diagram contributing to a and then "subtracting" all the physics contributing to the relevant lower 
scales. In the effective field theory (EFT) methodology this is equivalent to a multiple matching 
procedure starting from full QCD and the appropriate effective theories. In both approaches it is 
obvious that the hard part should depend only on the hard scale and the renormalization one as well. 
Thus, H is a polynomial of only one quantity: \og(Q 2 /fi 2 ). Any dependence on nUV regulators 
needed to calculate a is prohibited. 

2. The un-canceled nUV divergences in a, no matter how they are regulated, are completely a 
genuine IR ones that has to be generated by similar ones from the right-hand side of Eq. ( |1.1[ ). This 
is the ultimate check (at least pertubatively) for either the subtraction philosophy or the EFT match- 
ing procedure to work! The immediate conclusion is that the product J n JnS has to include only 
the IR of QCD. Again we emphasize that this simple observation is (or should be, if perturbative 
calculations are performed properly) independent of any regulator(s) of the nUV divergences. 
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Given the factorization theorem of a, where a is calculated in full QCD, with the above two 
observations and combined with the requirement of identical treatment among the two collinear 
sectors n and n, one can easily conclude that each quantity, J n y/S and J„ y/S, has to be a well- 
defined quantity in the sense that it is free of any nUV divergences that do not exist in the full 
pQCD calculation of a (a calculation that is needed in order to extract the hard part H). Again 
we mention that this line of reasoning is valid no matter which set of regulators is used for the 
nUV as long as one uses the same set for both collinear sectors. Thus we define the quantity 
F n (n) = Jn,(n) as the "TMDPDF". For two identical partons in the initial states one simply has: 

Notice that no attempt to go "off-the-light-cone" has been made. There are mainly two reasons 
for that. First is that one is completely free to regularize the nUV divergences with any set of 
regulators and moreover there is also no need to even distinguish between the different kinds of 
divergences in that set. This is definitely true if one believes that the factorized result of a holds to 
all orders in perturbation theory. The second reason is that large logarithms can be resummed (and 
one should be able to resum them) without any reference to a specific set of regulators. We believe 
this observation is simple enough and it should be undisputed. 

We finally remark that if one insists on introducing different regulators for the two collinear 
directions (like different off-shellnesses, A's, etc.), as it is done in [|3|], then certain logarithms of 
the ratio of those different « ^-regulators will still appear in each one of the TMDPDFs, F n and 
F w . Those logarithms are intimately related to the existence of the Wilson lines in the collinear and 
soft matrix elements. However those logarithms will cancel when one combines the two collinear 
contributions and the soft function to a. Thus we find it completely harmless to use identical 
parameters as long as one can still resum large logarithms and extract each TMDPDF from the 
relevant experimental data. 

We start by briefly reviewing Collins definition of the TMDPDF [Q] and then compare it with 
the "on-the-light-cone" one Some phenomenological applications of the definition and evolu- 
tion of TMDs as defined in ^ can be found in Ref. [0]. 



2. Collins Approach 

In impact parameter space, Collins defines the TMDPDF as: 



on-the-LC off-the-LC 

. / 

X 2 S(h'-i-oo v„\ 

F n (x,b;Q,y n ) =J n (x,b;Q ) 



S(b;+°°,y n ) UxMQ 2 ) S{b;+™,y n ) 



S(b;+oo 7 -oo)S(b;y n ,-oo) y/§(b; +00, -00) y S(b;y n ,-°°) 

(2.1) 

The first factor in the last term is exactly the definition of the TMDPDF where all Wilson lines 
are still on-the-light-cone. However it should be noted that this definition relies on the hypothesis 
that soft contamination in the naively calculated collinear contribution J n can be removed by soft 
function subtraction. This issue is a subtle one and its validity depends on the consistent use of 
the nUV regulators of the naive collinear, its soft limit (or its zero-bin) and the soft function built 
from soft Wilson lines. More discussion on this can be found in [f^, || [J. Hereafter we will not 
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discuss this issue further. The aim of the second term, is not canceling the rapidity divergencies 
as one might naively think (especially when going off-the-light-cone) but to introduce a new pa- 
rameter y n that allows the derivation of a "Collins-Soper" like evolution equation to resum large 
logarithms of Qjkj. The dependence (of Collins' TMDPDF) on this arbitrary parameter will can- 
cel only when two TMDPDFs (or, more generally, two "collinear" hadronic contributions like one 
TMDPDF and one fragmentation function in the case of SIDIS) are combined since the two sectors 
are tilted asymmetrically (notice the e y " in the two tilted collinear directions) so as to achieve this 
cancelation. We want now to raise two important points: 



1. The second factor in the last term of Eq. (2.1 ) introduces an arbitrary dependence on y„ of a 
single TMDPDF. In all existing literature using Collins formalism (e.g., [ 12, [13], |14], |15|]) for 
single sector, this dependence is canceled when y n is set to zero, which is to say that we are 
back on-the-light-cone since that factor reduces to unity. 

2. If that factor is exactly 1 on-the-light-cone, and one can still resum large logarithms by 
simply noticing that the first factor in the last term of Eq. (2T) is free from all un-physical 
nUV divergences (under the requirement of identical treatment of two collinear sectors) then 
why one needs to consider the y„ -dependent factor in the first place? Life is much simpler 
without it! 

Given all the considerations above, we now show more explicitly how Collins approach re- 
duces to ours when one chooses y„ = 0. In that approach the master equation for quark TMDPDF, 
denoted by F^, p , is 

/l A V f Y \ 

X 

where the matching of the quark TMDPDF onto the standard PDF is 

/l A-J x 
— C c f/j (-,b^F^)t/p(. x '^) + ^(( A QCDb) a ) . (2.2) 

J=*# X X 

On the other hand, our master equation for the TMDPDF on-the-light-cone is 

- 1 dx' 



j=q,g J * x Vx ' 

{ , , Of f^= Q fdu' 
xexp{-D{b;ii I )ln^+ 



where the matching of the quark TMDPDF onto the standard PDF is given by 

F f/P (x,b;Q 2 f ,n)= £ / — C f/j (-,b;Q 2 f^)<Pj/p(x^) + ^((A QCD br). (2.3) 

;— ^ „Jx X X 



j=q,g 
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The fact that the TMDPDF Fyyp is free from all un-physical nUV divergencies allows us to write 

c flj {xMQ\n) = CfjfcM. (2-4) 

In both cases should be chosen so that the large logarithms in the matching coefficients C 



are canceled. This choice can follow the well-known CSS approach fllOQ, where one relates /I with 
b and has to deal with the issue of Landau pole. However one can also follow the work of [ 11 ] and 
transform back to momentum space before setting the scale which is related finally to kj. 
When setting y n = one gets the following relations: 

K(b;n)^-2D(b;n), — ► y F {Q 2 ^) , Cfj/fab; Cm) — > C f/j (x,b;Q 2 ^) ■ 

(2.5) 

The above relations hold to all orders in perturbation theory. Thus we see that Collins formalism 
reduces to the "on-the-light-cone" one. 



3. The TMDPDF On-The-Light-Cone 

For many years, the two main issues that people have tried to resolve in order to get the correct 
definition of the TMDPDF are the following: 

• On one hand, the TMDPDF should be free from rapidity divergencies. In the case of the 
PDF, as it is shown in Sec. |], these divergencies cancel when one combines virtual and real 
diagrams. But when we allow the collinear matrix element to depend also on the transverse 
coordinate, such cancelation does not hold anymore. 

• On the other hand, the TMDPDF is a quantity that depends on two different scales: the hard 
probe Q and the transverse momentum kj. When these scales are widely separated, one 
needs to resum the large logarithms appearing in perturbative expression of the TMDPDF. 

In the following we give a consistent definition of the TMDPDF, while staying on-the-light- 
cone that addresses the issues mentioned above. More details can be found in [Bp. 

An intermediate step towards getting the factorization theorem of the hadronic tensor M is 
given by 

M = H{Q 2 /pl 2 ) j d^e-^y J n {0 + ,y- ,y i _)Hy + ,0- ,y i _)S(0 + ,0- ,y L ) , (3.1) 

where H is the hard matching coefficient, J n r n \ are the pure-collinear matrix elements and S is the 
soft function. The collinear pieces J n r n ^ are intended to contain just collinear modes, however in 
actual calculation we integrate over all momentum space thus collecting a contribution from soft 
modes as well. Below we denote the "naive collinear" as: J n rj(\. 

The fact that the soft function depends just on the transverse coordinate was derived in [Q] 
based on power counting arguments of the soft and collinear modes. 
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Figure 1: Virtual corrections for the collinear matrix element. The black blobs represent the collinear 
Wilson lines W in Feynman gauge or the T Wilson lines in light-cone gauge. Curly propagators with a line 
stand for collinear gluons. "h.c." stands for Hermitian conjugate. 



Based on Eq. (fO|), by symmetry, we define the TMDPDF as 

F n (x;k nX ) = y ^_^ e -(l^ + -a4x) / , i( o+ )r - ) ^)^ 5 (o+,0-,r x ) 

dr~d 2 r ± i{ l 2r - xp+ _ y , r kn , ) f» (0+ , r- , r ± ) 
(2*) 3 ^(O+.O-.rx)' ^ 

where the second equality holds only if the subtraction of soft contamination is equivalent to divid- 
ing by the soft function. 

We can use any regulator to regularize the nUV, however one has to be consistent. A factor- 
ization theorem can be understood from the effective field theory point of view as a multiple steps 
of matching between different effective theories. In this sense it is a must that one regularizes the 
"IR" consistently in the theory above and below the relevant scale. Then, the matching coefficients 
(H and C in our case) can never depend on the set of the nUV regulators. 

We choose a frame where p = (p + , 0~ , 0^ ) and p = (0 + ,p~,0±) with p + =p~ = Q, and write 
the poles of fermion propagators with a real and positive parameter A: 

+ # , + # + # , + # (33) 



(p + k) 2 + i0 (p + k) 2 + iA' (p + k) 2 + i0 (p + k) 2 + iA' 

The above prescription applies as well to the fermion propagators. The corresponding pole-shifting 
for collinear and soft Wilson lines goes as follows 

(3.4) 



k+±i0 k+±i8' k-±i0 k-±i5' 
where 8 is related to A through the large components of the collinear fields, 

« = 4 = =. (3-5) 

P+ p 

This relation comes from the fact that in order to recover the correct IR structure of QCD, the 
contributions of collinear and soft Wilson lines have to be consistent with the collinear and soft 
limits, respectively, of full QCD. 
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+ h.c. 



(a) 



(b) 



(c) 



Figure 2: Virtual corrections for the soft function. Double lines represent the soft Wilson lines, S n ^y "h.c' 
stands for Hermitian conjugate. 



3.1 Virtual Diagrams 



The diagrams in figs. ([I]) and (|2|) give collinear and soft virtual contributions respectively to 
F„. The Wave Function Renormalization (WFR) diagram (|I|a) and its Hermitian conjugate give 



(l-x)8^(k nX ) 



1 , At 2 1 
hln— + - 

£uv A 2 



(3.6) 



The W Wilson line tadpole diagram, ([]})), is identically 0, since n 2 = 0. Diagram (JT|c) and its 
Hermitian conjugate give 

7 k 2 ' 



In 



■In" 



£uv 



8A 



2hw+ln z ^-+2- 



12 
(3.7) 



The contribution of diagrams (|2|a) and (|2|b) is zero, since (|2]a) is proportional to n 2 = and (Qb) to 
n 2 = 0. The diagram (||c) and its Hermitian conjugate give 



a s C F 
2n 



8® (I, 



-ln^r+ln 2 ^r + 



7T 



-uv 



£uv M 



(3.8) 



The virtual part of the TMDPDF at 0{a s ) while using the relation in Eq. ( |3.5[ ) is 



r nl — 2 nl nl 2 



a s Cp 
2n 



8{\-x)8 {2 KKl) 



1 1 

— + ■ 

-uv 



3 , /I' 
+ ln" 



£uv V 2 G 2 / 2 M 



3, A 1, j A 2 2 A 7 7T 2 " 1 
: ln— - -In 2 _„ „ + ln 2 — H 



2,,2 



e 2 M 



(3.9) 



As mentioned earlier, individual contributions to F v . ,_. , have mixed divergences, however F\ itself 
is free from them. The A-dependence that remains is pure IR, as can be seen from the matching 
with full QCD in Sec. |. 



3.2 Real Diagrams 

The relevant diagrams for the real part of F n are shown in figs. (|J) and @. Diagram ([|a) gives 



/ f) = ^( 1 _ e)(1 



k 2 



\k 2 nT -/A(l -x) 



(3.10) 
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(a) 




(b) 



+ 




(') 



(d) 



Figure 3: Real gluon contributions for J 



n n 
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(a) 



(b) 



n n 
i 




(c) 



(d) 



Figure 4: Real gluon contributions for the Soft function. 



The sum of diagram (||b) and its Hermitian conjugate ^p) is 



'ni 



2n 2 



{\-x) + i8/ P + 



k 2 nT -iA(l-x) 



+ h.c. , 



(3.11) 



Diagram (||d) is zero, since it is proportional to n 2 = 0. 

For the real emission of soft gluons, diagrams (0a) and (Qd) are zero, since they are propor- 
tional to n 2 = and n 2 = respectively. Diagram (|]b) and its Hermitian conjugate ^p) give 



a s C F 



% 2 k 2 nT -8 2 "'k 



i 52 
In—. 



(3.12) 



Combining the above contributions according to Eq. ( J3.2| ), one gets the real part of the TMD- 
PDF F^j at 0{tts), which is free from rapidity divergencies. A detailed discussion about the re- 
maining pure IR A-dependence can be found in [S]. 



4. Integrated PDF 

In this section we briefly report the calculation of the integrated PDF at first order in a s with 
the A-regulator. As it is shown below, the mixed divergencies will cancel when we combine virtual 
and real contributions, leading to the well-known IR collinear divergence of the PDF. 

The virtual diagrams for the PDF are the same as for the naive collinear matrix element that 
enters into the definition of the TMDPDF, Fig. (fil). Then we have 



^1 = ^5(1 



1 1 /3 , H 2 \ 3, A 1 2 A 2 2 A 7 n 2 ^ 
■0- + ;7 + ln^ -^ln^---ln i — — + ] n 1 — + -- — 



-uv 



£uv V2 ' "'Q 2 J 2 ]X 2 2 Q 2 ]X 2 \x 2 ' 4 3 



(4.1) 
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The real diagrams are the same as in Fig. (|3|), from which we get 

.2 



2n 



— + ln— -l-ln(l-jc) 
£uv A 



(4.2) 



and 



, (||+|) = (XsCf 
Vnl 2% 



( 1 At 2 
+ ln^- 

V £uv A 



<rV 26(1 -*l 



-25(1 -x) 1 



7T 2 1 



-In 2 



7T 



24 2- fi 2y+Y 5 ( 1 -^) 



(4.3) 



Diagram (||d) is zero since it is proportional to n 2 . 

Mixed divergencies ^Mn (and also double IR poles ln 2 ^r) are cancelled when we combine 
virtual and real contributions. The PDF to first order in a s is finally 



(j> nl (x;n) = 8(l-x) + 



a s C F 

271 



1 



1 



" qlq Vi^ _1 V ) ~ 4 5(1 ~ X) ~ (1 ~ X) [1 +ln(1 ~* )] 



(4.4) 



where is the one-loop quark splitting function of a quark in a quark, 



1+x 2 
l-x 



l+x2 ■ 3 £M ^ 
-oil —x 



(l-x). 



7 -*+(l-x) + U{l-x). (4-5) 

^1— XJ + L 



The IR collinear divergence is encoded in the single logarithm In ( A/jU 2 ) , while the rest is a function 
of x that in general depends on the regulator used. 

5. Q 2 -Dependence and Resummation 



When kj S> Aqcd we can perform an operator product expansion (OPE) of the TMDPDF onto 
the integrated PDF, 



r 1 dx' / x \ 
F n (x;b,Q,n)= / —C n (-;b,Q,n)<t) n (x;n} 



(5.1) 



where 



i / , i Mr 

UX ^ ) = 2J^ 



. e -ib-*P + {p\Xn(0 + ,y-,0AxZ(0 + ,0-A)\p)U included- 



(5.2) 



and C„ is the matching coefficient that cannot depend on any IR regulator as mentioned before. In 
this section we compute C„ to first order in a s . 



The virtual part of the TMDPDF in momentum space was given in Eq. (4.1), and in impact 
parameter space it reads 



2% 



1 



+ • 



1 



uv 



- + In I — ^ln-^r — -!-ln^ 



21 



£uv V2 ' ~Q 2 ) 2"> 2 T" Q 2 ]X 2 



, 7 A 7 71 

+ln V + 4"T 



(5.3) 
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The Fourier transform of real-gluon emission diagrams given in Eqs. (3.10, 3.11, 3.12), while 
keeping the A's to regulate the nUV divergences, are 



1 — x)ln- 



2k 



A{l-x)b 2 ' 



(5.4) 



A 



a s C F 
2k 



ln- 



4e -2r E 



2.x 



Ab 2 \(l-x)- 

. 2 d(l-x)(\- — --ln 2 ^)} 
v ' V 24 2 Q 2 )\ ' 



A 

"0 



25(l-x)ln^2 ) + ^-5(1 - v) 



T 



and 



, _ a s C F f^ 2 4e- 2 KQ 2 | 2% 2 



2k \"' A 2 b 2 3 
Finally, the TMDPDF in impact parameter space to first order in a s is 



(5.5) 



(5.6) 



F n = <Pn + 



a s C F 
2k 



Lr-A, „ + < 1 -v) 5(1 -x) [ X -L\ - Ilj+LtI^ + ] o 



(5.7) 



Q 2 , 7T 2 
1— 

where Ly = ln(/x 2 fe 2 e 2r£ /4), 0„ is the PDF given in Eq. (4.4) and the remaining part is exactly the 
OPE matching coefficient 

i2 -2 \ 1 

7" T. J — 

1 2 

(5.8) 



C„(x;fc,e )J u) = S(l- JC ) + 



a. s C> 

27T 



-^L T + (1 -x) - 5(1 -x) ( X -L\ - ^L T + ln^L T + *L 



In the integrated PDF the mixed divergencies are canceled when we combine virtual and real 
diagrams. When we allow for a transverse separation also in the collinear matrix element, real 
diagrams do not cancel these divergencies of the virtual part. It is the soft function (square root) 
that does this job. Mixed divergencies in the virtual part of the collinear matrix element are canceled 
by the virtual part of the soft function; and the ones in the real part of the collinear are canceled by 
the real part of the soft function. All the remaining A-dependence is pure IR. 

The fact that the TMDPDF is free from rapidity divergencies allows us to extract and expo- 
nentiate its £2 2 -dependence to all orders in the following way [|5j], 



F n (x;b ± ,Q,n) 



Q2 b 2 \ -D(a s ,L T ) 



4 e -2y E 



C%(x;bx,n)®<i>n(x;n), 



(5.9) 



where 



C^(x;b x ,n) = 8(l-x) + 



a s Cf 
2k 



q/q L T + (1 -x) - 5(1 -x) ( --L 2 -^ T + 7 L 



(5.10) 



Given the renormalization group invariance of the hadronic tensor M in impact parameter 
space, 

M = H(Q 2 /ll 2 )F n (x;b ± ,Q,ii)Fjr(z;b ± ,Q,ii) , (5.11) 
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we can establish the following relation between the anomalous dimension (AD) of the hard match- 
ing coefficient, and the one of the TMDPDF, y„, 



Yn 



1 



1 



Q 



l 



~Yh = --A(a. v )hw - -B(a s ) 



(5.12) 



Notice that this relation allows us to extract the AD of the TMDPDF up to third order in a s from 



the known Jh at three loops [Q |19[ ]. And since A(a s ) = 2T cusp (a s ) to all orders in perturbation 
theory, we get also 



dD(a s ,L T ) 
dln/j, 



cusp 



a s ). 



(5.13) 



Our master equation for the quark unpolarized TMDPDF Ffj P m impact parameter space is 

- 1 dx' - 



F f/P (x,b;Q 2 ,ll = Q) = [ —Ctjf^M^Qj/ptfWi) 



J=9,g 



Q 2 fV=Q d\x' 



xexp\-D(b,^)ln-^+ — 



Q2 



(5.14) 



6. Hard Part at 



.a. 



The hard matching coefficient for the qj -dependent DY cross section is the same as the one for 
inclusive DY. As mentioned before, this matching coefficient at the higher scale Q is obtained by 
matching the full QCD cross section onto the imaginary part of the product of two effective theory 
currents. This echoes the "subtraction method" in perturbative QCD. 

We start by rewriting the cross section in a more useful way, 

4na dxdzd 2 q±^ 2 „w - ^ 



M(x,z:,q±,Q)=H(Q 2 /ll 2 ) J d 2 k n± d 2 k^8 {2 \q ± -k n± -k Tll _) [8(\-x)8^(k n± )8(\-z)8^(k^) 
+a s (F n{ 8(l-z)8^(knx) +F m 5(1 -x)8^(k nl _)^] +0(a 2 ) 
= H(Q 2 /ii 2 ) [8(l-x)5(l-z)8^(q ± ) 

+a s (8(l-z)F nl (x;q ± ,Q,n) + 8{l-x)F m (z;q ± ,Q,n))] +0{a 2 ), (6.1) 



where M is the hadronic tensor. 

In QCD the virtual part of M with the A-regulator is 



M, 



QCD 



a s C F 
In 



8(\-x)8(\-z)8 {2 \q A _) 



21 



21n 2 -^ - 31n— r - - + — 



Q 1 



(6.2) 



The above result can be simply obtained by considering the one-loop correction to the vertex dia- 
gram for — )• Y ', with the inclusion of the WFR diagram while using the fermion propagators in 
Eq. (ph. 
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The virtual part of F„ at one-loop is given in Eq. (4. 1 ), and we have an analogous result for F%. 



Using Eq. (|6JJ) the total virtual part of the hadronic tensor M in the effective theory is 



Mscet 



2k 



)S(1 



-21n 2 



1 



+ ■ 



uv 



1 



3 + 21n 



Q 2 



Q 2 



A IX IX 

3ln Q 2+3ln! t 2+ Q 2 + 2 



2k 2 



(6.3) 



where the UV divergences are canceled by the standard renormalization process. We notice that the 

are the same, as they should, thus the matching coefficient 



IR contributions in Eqs. ( |6.2[ ) and 
between QCD and the effective theory at scale Q is: 



H(Q 2 /^ 2 



1 + 



a s C F 
2k 



-31n 



e 2 



In 



2E 
Q 2 



: + ■ 



lK 2 



(6.4) 



The above result was first derived in [2C]. We can also obtain the AD of the hard matching coeffi- 
cient at ff{a s ) and verify Eq. ( |5.12[ ), 



oc s Cp 
2k 



-2y n i 



(6.5) 



It is clear that the IR divergences of full QCD are recovered in the effective theory calculation, 



Eq. ( p.3| ). To get this one must use the same regulator in both theories and the way that one 
regulates the IR physics (or more generally, the nUV) in both sides have to be consistent. The 
matching coefficient at the higher scale depends only on the hard scale Q 2 as it should be. 



7. Gauge Invariance: T-Wilson Line 

One of the advantages of setting all the perturbative calculation on-the-light cone is that it is 
straightforward to show the gauge invariance of the relevant physical quantities. Below we derive 
the transverse gauge link, the T Wilson line, that has to be included in the effective theory in order to 



accomplish the calculation of collinear and soft matrix elements in a gauge invariant way [21 



First we recall some of the features of the gluon fields in QCD in light-cone gauge (LCG) [p4|]. 
To fix matters, we work in QCD with the gauge fixing condition nA = 0. The canonical quantization 
of the gluon field proceeds by inserting in the Lagrangian the gauge fixing term J£ g f = A a (nA a ). 
The A" is a field whose value on the Hilbert space of physical states is equal to zero. It is possible 
to write the most general solution of the equation of motion of the boson field A^ by decomposing 
it as 



= T^k)8(k 2 ) +n ll ^-A%nk,k ± ) + l -^8(nk)U a (nk,k ± ) 



(7.1) 



where the field Tfi is such that n^T^ik) = and k^T^(k) = 0. Fourier transforming this expression 
we see that in general the field A^(jc) has non- vanishing "— " and "_L" (respectively nA a (x) and 
Ay a (x)) components when x~ — > ±°°. Now we define 



def 



(x ,xx) = A(x + ,oo , x± ) 

r,x±) 



(7.2) 



A(x + ,x ,x±) d = A(x + 
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and get 



np L = q ± + 8 4 _l = id/ L + ^ + g 4^ = f iij> x + g4^ . (7.3) 



Taking into account that the fields A± and in Eq. ( fL4| ) are evaluated at space-like separated 
points then one can show that 

Up± = Tify ± T\ (7.4) 



where 

T 1 = Pexp 



-ig [ dzl± •A ( L oo) (x+,x_ L -l±r) 
Jo 



(7.5) 



This equation leads us automatically to include the r -Wilson line at the level of the soft- 

collinear effective theory (SCET) Lagrangian. Moreover l/nd and T commute because the T- 

iH. 



Wilson line does not depend on x . Under gauge transformation 8A^j_ = D^^co one has T (x + ,x± ) 

iH 



U(x + ,xj_)T(x + ,xj_)U^ (x + ,xj_ — 1±°°) = U(x + ,xj_)T(x + ,xj_) since A^°°)(x + ,°°x) =0. Notice also 



that the r -Wilson lines are independent of l±. 

Now we split the fermion field into large and small components using the usual projectors 
and rjtijt/4 and eliminate the small components using the equations of motion [25]. The result 
of this is 

JSf = l n (inD + ilp^Ip^j . (7.6) 
In QCD and in LCG with the gauge condition nA = we get 

JSf = l n (inD + Tiip^iip^ . (7.7) 
In order to get the SCET Lagrangian we must implement multipole expansion and power counting 



on the fields that appear in Eq. (7.6). In SCET we have also the freedom to choose a different gauge 
in the different sectors of the theory. We distinguish the cases of SCET-I and SCET-II. The two 
formulations differ essentially in the scaling of the soft sector of the theory. Here we are interested 
in SCET-II which is the effective theory necessary for the Drell-Yan qj -spectrum where qj is much 
larger than Aqcd and much smaller than Q . The collinear modes we have in In SCET-II scale as 
(nk,nk,kx) ~ 2(1,T] 2 ,T]) where t] <^ I and the soft modes as (nk,nk,k±) ~ Q{r},r},r}). Collinear 
and soft gluon fields' components scale accordingly. 

The soft sector of SCET-II has peculiar features. In regular gauges soft partons do not in- 
teract with collinear ones because the interactions knock the collinear fields far off-shell. This 
is also true in light-cone gauge except when one makes the choice nA s = (take here a covari- 
ant gauge for collinear fields for fixing ideas). It is easy to be convinced that interactions like 
Yli<Pn{x)A'^_ L (x ~ ,x L ), where here "oo" refers to the + direction and (j> n (x) are generic collinear 
fields, preserve the on-shellness of the collinear particles. Using multipole expansion the ver- 
tex becomes Yli<l> n (x)A^(0^ ) (because for collinear fields x~ ~ 1 and for the soft field x~ ~ 
1/t]). In this gauge the covariant derivative for collinear particles becomes iD^ = id* 1 +gA„(x) + 
gA^^(0~ ,x±). The gauge ghost A^ however can be decoupled from collinear gluons defining a 
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"soft free" collinear gluon A^ 11 (x) = T sn (x±)A% (x)T^ n (x±) where 



T m =Pexp 



dzl±-A s ^\o ,x± — 1±t) 



(7.8) 



Defining 



id* +gA ( n ]tl we have il ± D ± = T n {x±)il±D^T^(xx) and 



(7.9) 



where = and W n T = rJ ^W^ ' are made out of soft free gluons. Thus, and thanks 

to T sn Wilson lines, the soft partons are completely decoupled from the collinear ones. 

With the existence of the r -Wilson lines it is possible now to compare the one-loop contri- 
butions in both Feynman gauge and LCG. This comparison can be performed at the level of the 
integrands without the need to perform any actual calculations. 

In light-cone gauge we use the ML prescription [23 ], which is the only one consistent with the 
canonical quantization of QCD in this gauge [24]. Moreover in the n and n collinear sectors the 
only gauge fixings compatible with the power counting of the collinear particles are respectively 
nA n = and nAn = 0, which correspond to "killing" the highly oscillating component of the gluon 
field in each sector. We now compare the integrals that we have evaluated in Feynman gauge with 
the corresponding ones in light-cone gauge. 

The interesting contribution to the collinear part of the TMDPDF in Feynman gauge is pro- 
vided by the W n Wilson line and it is (cfr. Eq. (3.7)) 



^) (Feyn) 



8{l-x)8^{k n± )2ig 2 C F n 



2f 



d d k 



1 



p++k+ 



(2n) d (k 2 + iO) (k+ - iO) (p + k) 2 + iO ' 
In light-cone gauge this result is reproduced when one combines the axial part of the WFR 



(7.10) 



ffrHM = 8{\-x)8^Ck n xWc F ^ 



1 



p++Jt H 



d d k 



{2n) d (k 2 + i0) (p + k) 2 + i0 



+ 



d(-k-) 



e(k~ 

k+ + i0 ' k+-i0 

(7.11) 



with the contribution of the T Wilson line 



=-8{\-x)8^CKx)2ig 2 C F ^ 



d d k 1 p + + k + 
{2n) d {k 2 + i0) (p + k) 2 + iO 



i 



iO k+ + i0i 
(7.12) 



It is clear that f§' ] {Feyn) = f§ :) {T) - f$ ,] {Ax) /2. In SCET there is also a tadpole diagram that 
needs to be considered however it is null in light cone gauge since the gluon field does not propagate 
at infinity [21]. Similar considerations hold for Feynman diagrams with real gluon contributions 
and for the soft sector [j|]. 

As a final remark, it is not so obvious how one can implement light-cone gauge calculations 
while going off-the-light-cone. 



14 



TMDs 



Ignazio Scimemi 



8. Conclusions 

We have argued that to properly define the TMDPDF, one needs to include the soft factor (a 
square root of it). With this inclusion the TMDPDF can be renormalized and it captures all the IR 
of full QCD when it is introduced in a properly obtained factorization theorem. We argued that one 
does not need to introduce any of the "off-the-light-cone" contributions, while still being able to 
resum large logarithms, thus avoiding the use of the Collins-Soper evolution equation. 

The evolution of the TMDPDF has been discussed and the <2 2 -resummation at the intermediate 
scale was explained. That was achieved under the assumption of identical treatment of two collinear 
sectors and it was mainly based on the fact that the TMDPDF is free from rapidity divergences. 

The conclusions above are independent of the set of regulators that one uses for the nUV. 
Moreover, if one is consistent with the regularization of nUV divergencies in the theories above 
and below a certain relevant scale, then the matching coefficients will be independent of this set of 
regulators. 

We also discussed the relevance of the r -Wilson line and their contribution to render hadronic 
matrix elements gauge invariant. Contrary to many existing statements in the literature, those 
Wilson lines are crucial for the gauge invariance to hold. Finally, all this effort can be extended 
to the proper definition and treatment of polarized quark and gluon TMDs, which are important 
ingredients in order to get information about the inner structure of the nucleons. 
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